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ABSTRACT 

We present new measures of the evolving scaling relations between stellar mass, luminosity and 
rotational velocity for a morphologically-inclusive sample of 129 disk-like galaxies with zab < 22.5 in 
the redshift range 0.2 < z < 1.3, based on spectra from DEIMOS on the Keck II telescope, multi- 
color HST ACS photometry, and ground-based Ks-band imaging. A unique feature of our survey 
is the extended spectroscopic integration times, leading to significant improvements in determining 
characteristic rotational velocities for each galaxy. Rotation curves are reliably traced to the radius 
where they begin to flatten for ~90% of our sample, and we model the HST-resolved bulge and 
disk components in order to accurately de-project our measured velocities, accounting for seeing and 
dispersion. We demonstrate the merit of these advances by recovering an intrinsic scatter on the stellar 
mass Tully-Fisher relation a factor of 2-3 less than in previous studies at intermediate redshift and 
comparable to that of locally-determined relations. With our increased precision, we find the relation 
is well-established by (z) ~ 1, with no significant evolution to (z) ~ 0.3, AM* ~ 0.04 ± 0.07 dex. A 
clearer trend of evolution is seen in the B-band Tully-Fisher relation corresponding to a decline in 
luminosity of AMb ~ 0.85 ±0.28 magnitudes at fixed velocity over the same redshift range, reflecting 
the changes in star formation over this period. As an illustration of the opportunities possible when 
gas masses are available for a sample such as ours, we show how our dynamical and stellar mass data 
can be used to evaluate the likely contributions of baryons and dark matter to the assembly history 
of spiral galaxies. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: kinematics 
and dynamics — galaxies: spiral 



1. INTRODUCTION 

A major challenge for ACDM structure formation lies 
in understanding how the baryonic components of galax- 
ies assemble within dark matter halos. Although baryons 
represent only one-s ixth of the gravitat i ng matter in 



WMA P cosmology ( |Spergel et al.||2007| Seljak et al 
2005 1 , their dissipative properties sugg est they dominate 
the in ner regions of luminous galaxies (Blumenthal et al. 
19861. Determining the interplay between dark matter 
and baryons is critical for predicting the evolution of 
density profiles, subs tructure, shapes, and angular mo 
mentum of galaxies ( |Governato et aL 2007 Shlosman 
2009|) . One of the most significant challenges is reproduc 



ing the detailed characteristics of rotationally supported 
disk galaxies which represent the dominant frac tion of 
present-day luminous systems ( Ellis fe Silk||2009 1 . 



Observational efforts in this challenge have focused on 
the Tully-Fisher (TF) relation fTufly fc Fisher] [1977 1 
and its past evolution. This important scaling relation, 
which correlates disk luminosity with rotational velocity, 
provides an essential benchmark for verifying theoretical 
models based on the standard dark matter picture. Early 
N-body simulations as well as semi-analytic models pro- 



van den Bosch 


2000 Mo & Mao 


2000 Eke et al. 


2001 



s .miller 1 ©physic s. ox. ac.uk| 

1 Oxford Astrophysics, Oxford, OX1 3RH, UK 

2 Astronomy Department, University of California, Berkeley, 
CA 94720 

3 California Institute of Technology, Pasadena, CA 91125 

4 UC Santa Barbara Physics, Santa Barbara, CA 93106 



Benson et al. 2003 Dutton et al. 2007). Caused by a 
transfer ot angular momentum from baryons to the dark 
matter halo, this deficiency has since been mitigated by 
impro ved resolution, as well as th e introduction of feed- 
back (jSteinmetz & Navarro 199 9]), e.g. from supe rnovae 
( |Governato et al.|2007[ JPiontek fc Steinmetz|2009[ ) . How- 
ever reproducing the absolute values observed in the scal- 
ing relation has remained problematic. 

Despite these challenges, the theoretical understanding 
of disk galaxy scaling relations and their evolution has 
made some improvement over the past decades. Using an 
adjustment to the rotational velocity derived from their 
hydrodynamic simulations to accou nt for over-merging, 



Portinari & Sommer-Larsen (2007) were able to match 
the observed local TF relation, and claim a modest evo- 
lution to z 1. However the predictive power is tem- 
pered by an unknown dependence on redshift of this 
adjustment. Semi-analytic models have also worked to 
match observations and provide further insight on the 
physical interpretation of evolution in the TF relation. 
Some controversy remains over whether the central re- 
gions of galaxy halos ar e subject to adiabatic contraction 
(Somerville et al. 2008), broadly ma intain a non-evolving 
density pro file ( Wcchsle ret al. | 2 002), or permit adiabatic 

"TBI 



expansion flDutton et al.|2010b[ ). Regardless of the exact 



evolutionary response ot the inner halo, the persistent 



in tandem with the dynamical mass dFall & Efstathiou 


1980 


Dalcanton et al. 


1997 


Mo et al. 


1998 


). Gas may 



increasing the disk scale length as stars form. In this 



2 



Miller et al. 



framework, while any given galaxy is predicted to grow 
by factors of 1.2-2 in stellar mass, dynamical mass, scale 
radius, and luminosity since z ~ 1 (modulo evolutionary 
corrections), this growth typically occurs along scaling 
relations, reducing the evolutionary signals accessible to 
observations. 

Observational progress in testing these pictures of disk 
assembly has been similarly slow. There are significant 
technical challenges in making the necessary measure- 
ments at intermediate redshift and, as a result, there 
are discrepant conclusions with regard to evolutionary 
trends in the literature. In part, this may reflect dif- 
ferent ways i n which intermediate red shift disk galaxies 
are selected. Vogt et al. (1996 19971 undertook an im- 
portant pioneering study, finding a modest increase in 
luminosity (AMg ~ 0.6) at fixed velocity to z ~ 1, but 
deduced this represented only an upper limit to possible 
evolution because of sample biases and other assump- 
tions. Subsequent optical-based studies have presented 
mixed conclusions. A key uncertainty is whether to ad- 
dress evolution in the overall mass-to-light ratio indepen- 
dent of lumi nosity (i.e. a zero - point shift with redsh ift) as 



discussed by Rix et al . . 
|Fernande z Lorenzo et a. 
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aOT2009. 
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and 

' PT097pnlO[ ), or whelherto per- 
mit luminosity dependen t evolution (i.e. cha nges in the 



TF slope) a s discussed by Ziegler et al. (2002) and Bohm 



et al. (2004). TF studies at infrared wavelengths are less 



affected by biases induced by short-term star formation 
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Clearly the rest wavelength a t which the luminosity is 
sampl ed is a key parameter: [Fernandez Lorenzo et al.| 
(2010) claim evo lution in the -B-band bu t none in red- 
der bands, while Weiner et al. I ( |2006a|b[ ) find evolution 
in the slope of the infrared TF relation consistent with 
that seen in the blue relation, h owever they observe little 
evolution in infrared zero-point. Moran et al. ( 2007[ ) have 
emphasized the importance of environmental influences 
which can be deduced by considering the scatter in the 
TF relation as a function of local density. 

In view of this, a more physically-relevant approach 
for understanding the assembly history of disks may be 
to consider the stellar mass TF relation (M*-TF) which, 
notwithstanding the difficulty of estimating gas fractions, 
provides the most robust route to understanding the in- 
terplay between baryons and dark matter in disk galax- 
ies. Stellar masses are derived using population model 
fits to multi-color photometry for galaxies of kn own red 



shift, assuming an initial stellar mass fu nction (Brinch 



mann fc Elhs||2000| [Bundy et al.||2005[). Altho ugh the 



low redshift M^-'L'b' relation is well-constrained (Bell fc 



de Jong| 2001| |Pizagno et al |[2005 j |Meyer et ~al. 
those at inter mediate redshift |(Jo^eliceeF^rpUDo| ' I 1 'lo- 



res et al.|2006 Atkinson et al.|2007 ) reveal a larger scat 



ter than seen in the traditional I F relations, suggestive of 
additional uncertainties. A recurrent topic of discussion 
in the literature is sample selection and whether evolu- 
tion seen in both the TF relation and its scatter is driven 
by redshift-dependent selection criteria. The inclusion 
of more early types and dynamically disturbed galaxies 
likely broadens the intrinsic scatter. In an attempt to 
include more kinematically-disturbed galaxies, often ex- 



cluded witho ut good cause in TF studies, | Weiner et al.| 
(2006b) and Kassin et al. ( 2007[ ) have included an ad- 
ditional kinematic term, 6 0.5, which accommodates the 
isotropic velocity width of the observed emission lines 
and reduces the scatter from their classic M*-rotational 



velocity relation. Kassin et al. (2007) detect no signifi- 
cant evolution in the IT 1 relation, including the So. 5 rela- 
tion, over 0.1 < z < 1.2. However, combining a measure 
of the velocity dispersion with the disk angular momen- 
tum may obscure information about the rotational sup- 
port of the disk. 



Using inte gral field unit (IFU) sp ectrographs, Flores 



et al. (2006) and Puech et al. (2008) have produced in- 
tensity, velocity, and dispersion maps of galaxies at inter- 
mediate redshifts that demonstrate the unique advantage 
of a second spatial dimension in modeling the velocity 
field and accounting for projection effects. So far, the 
IFU-based samples are fairly modest in size and sample 
brighter sources compared to those reached with multi- 
slit techniques. Moreoever, the spaxel resolution is often 
lower than for the highest-resolution slit spectroscopy. 
As we will show in this paper, the spectroscopic signal-to- 
noise is an equally important factor in making progress 
because it determines the radial extent to which emission 
lines can be traced and as a result, whether the adopted 
velocity measure requires extrapolation. 

Both IFUs and slit spectroscopy have been successfully 
employed at z sa 2—3, where near-IR studies can take ad- 
vantage of redshifted Ha lines — the brightest kinematic 
tracer — as well as improved spatial resolution from adap- 
tive optics. Results at these redshifts may indicate the 
emergence of regular scaling relations from more complex 
and disordered dynamical states. With adequate sam- 
pling, at least one-third of z ~ 2-3 star-forming galaxies 
show ordered rotation (Shapiro et al.|j2008| |Stark etaL 
|2008| I Jones et al. 2010 ) and reveal significantly higher 



velocity" disp ersions tha n loca l count erparts (e.g., Genzel 



et al. (2006)). Cresci et al. (2009) have measured the 
stellar mass TF relation at z ~ 2 using SINFONI ob- 
servatio ns of 18 rotation-dominated s ystems in the SINS 



survey ( Forster Schreiber et al.|2009[ ). The slope of their 
measured relation is consistent with that seen in local ob- 
servations but offset towards lower M* at fixed velocity 
by ^0.5 dex. While necessarily biased toward massive 
systems with well-ordered rotation, these may be repre- 
sentative of gas rich systems in transition to z^l disks 
( Tacconi et al.|2010 ). Gnerucci et al. (|2010[) have recently 
measured the TF relation at z ~ 3 from SINFONI IFU 
data, but because of the large scatter observed, they sug- 
gest the TF relation has not yet been established. How- 
ever all points on the relation are consistent with the 
favored vector of disk assembly theory, with a lower av- 
erage stellar-to-dynamical mass ratio than found in the 
local universe. 

The present survey was motivated by our desire to 
chart and understand this evolution from a prevalence 
of disturbed and complex dynamical states observed at 
high z to the well-ordered rotation of local spirals. To 
make progress, we seek to determine the stellar mass TF 
relation over the redshift range 0.2 < z < 1.3 with spec- 
troscopic exposures 3 to 8 times that of previous stud- 
ies. By including disk systems selected from HST ACS 
data with irregular or distorted morphologies, we hope 
to avoid biases based on selecting mature, well-ordered 
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disks ( Vogt et al.|19 96, 1997). Our study not only allows 
us to chart evolution m a large sample down to fainter 
limits and masses than is possible at z ~ 2, but the im- 
proved precision we demonstrate enables us to measure a 
robust TF relation only a few Gyr later. The scatter we 
observe should provide a valuable indication of the rate 
at which disks settle onto the local TF relation. To fully 
utilize the gains in signal-to-noise ratio from long expo- 
sures, we develop an improved method for extracting the 
rotation curves of galaxies at intermediate redshift via 
a new modeling code. While our results are based on 
slit spectroscopy, unlike most previous slit-based studies 



galaxies, i n the GOOD S Nort h and South fields, dis 
cussed by Bundy et al. (2005). A key diffe r ence with 



( Werner et al.|2006a| |Kassin et al.|2007| |Bohm fc Ziegler 
2UU7p , we are able to align the spectral slits on our masks 
with the HST-measured major axis, significantly improv- 
ing the fidelity of our recovered rotation curves at z ~ 1. 
We aim to avoid introducing a bias towards aligning the 
often more extended and brighter objects at lower red- 
shift over the higher redshift objects observed at smaller 
angular scale, since doing so could introduce an extrane- 
ous evolution in the offset and scatter of the Tully-Fisher 
relation with redshift. 

A plan of the paper follows. In §2, we describe our 
sample, the Keck spectroscopic data and the HST ACS 
resolved photometry in the north and south GOODS 
fields. Noting the limitations of earlier work, §3 intro- 
duces a new procedure for the analysis of rotation curve 
data. We justify our chosen method, discuss error esti- 
mation, and compare to previous work. In §4 we present 
the various TF relations, and in §5 we discuss methods 
for deriving dynamical masses to compare to baryonic 
mass estimates for a physical interpretation of our re- 
sults. Finally, §6 summarizes the overall conclusions for 
the assembly history of disk galaxies. 



Throughout the paper we adopt a Chabrier ( |2003 1 ini- 
tial mass function and a 57a = 0.Y, \l m = 0.3, Hq = 
70 km sec -1 Mpc -1 cosmology. All magnitudes refer to 
those in the AB system. 

2. DATA 

A prerequisite for constructing a disk galaxy sample 
suitable for measuring the evolution of the Tully-Fishcr 
relation is deep HST imaging essential for morphological 
selection, resolved photometry and accurate disk posi- 
tion angle data for the multi-slit spectroscopic campaign. 
Both north ern and southern GOODS fields (Dickinson 
et al. 2003[ ) are visible from the Keck observatory and 
remain the most appropriate areas for such a study given 
the unique availability of deep multi-color ACS data. In 
t his s ection we introduce our sample selection criteria 
(i ]2.ip and the spectroscopic data used for measuring the 
internal dynamics of our sample (j ^2.2[ ). We also intro- 
duce the photome tric data used for measuring galaxy 
sizes and sha pes (i 2.3 1 and stellar mass estimates from 
SED fitting (|2~4|. 



2.1. Sample Selection 

Our goal in the morphological selection of disk targets 
is to be inclusive of all galaxies with disk-like structure, 
avoiding the temptation of selecting the most relaxed and 
"well-behaved" spirals in favor of a more complete cen- 
sus, including the more disturbed and morphologically 
abnormal population. Sources were selected visually by 
co-author RSE from a zf$50LP < 22.5 sample of 2978 



earlier work (i.e. Vogt et al. ( 1996[ 19971; Conselice 



et al. ( 2005 )) is the inclusion of less mature morphological 



types which contain some evidence of disk-like structure 
as well as systems that may be interacting. We included 
visually irregular systems with elongated features, galax- 
ies with asymmetric and clumpy light distributions. We 
also included disks with dominant bulges. The main aim 
of broadening the morphological selection criteria was 
to avoid potential biases associated with selecting only 
symmetric spirals which may represent the end point of 
evolution and consequently bias us to locating mature 
systems. Within this z-band limited sample, we applied 
a further photometric selection, K s < 22.2, to ensure a 
high frac tion of reliable stellar mass (M*) estimates to 
(see {2A\, resulting in a morphologically suitable sam- 
ple of 1388 objects. Although spectroscopic redshifts are 
available for many of our targets from the Team Keck 
Redshift Survey program ( |Wirth et al.[2004| in GOODS- 
N and from the VIMOS VLT Deep Survey ( |Le Fevre 
et al.||2004| ) in GOODS-S, we did not exclude targets for 
which only photometric redshifts were available. Select- 
ing within our target redshift range of 0.2 < z < 1.3, we 



used phot ometric redshifts fro m COMB017 (Wolf et al . 
2004| and 



Bundy et al.| ( |2005[ ). As discussed later (j pL5| 
several galaxies within o ur sample can be f ound in ear- 
lier kinematic surveys of Flores et al. ( 2006 ) and Weiner 
et al. | (12006a I. 



2.2. Spectroscopic data 

Over a number of seasons we collected spectroscopic 
data for this sample with the DEIMOS (DEep I mag- 
ing Multi-Object Spectrograph; Faber et al. (2003)) in- 
strument on Keck II. In total, we examine d 23 6 galax- 
ies drawn from the target list discussed in [2.1 (17% of 
total sample), simply chosen to maximise number of ob- 
jects with best position angle alignment with that of the 
mask. Of the 236 spectra, 129 show line emission ex- 
tending past what we will term the "seeing-dispersion 
beam'12 59 have only very compact emission lines that 
sample the central region, and 48 are in a category we 
will refer to as "passive" (meaning galaxies which show 
no significant emission lines across the 2D spectrum, al- 
though weak lines may be recovered in the integrated 
spectrum) (Fig. [TJ Table [TJ. 

The bulk of our analysis is thus based on the 129 galax- 
ies with extended line emission. While this subsample 
makes up only 55% of our initially targeted sample, there 
is no statistical difference in its redshift or apparent mag- 
nitude distribution from the compact and passive subsets 
(Fig. [I]). The only obvious differences amongst the 3 
subsamples concerns their disk sizes and stellar masses. 
Disks with extended line emission have larger scale radii 
than those of the compact and passive subsamples, and 
passive galaxies are largely drawn from the upper end of 
the total stellar mass distribution (see Table [l]) . These 
differences are not surprising and do not lead us to sus- 
pect that our working sample of 129 disk galaxies is sig- 
nificantly biased in its range of physical properties com- 

5 The seeing-dispersion beam signifies a 2D Gaussian represent- 
ing the combination of the effect of seeing along the spatial axis 
and the emission line velocity dispersion along the wavelength axis. 
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Fig. 1. — The properties of galaxies in our sample in terms of the distributions of stellar mass estimates, spectroscopic redshifts, ZF850LP 
apparent magnitudes, and disk scale radii. Each histogram is partitioned according to galaxies which have extended emission lines (129), 
passive spectra with no emission (49) and spectrally compact sources where emission does not extend beyond the central-most regions of 
their disks (59) (see jj2~2|and Table [Tj for details). 



pared to the original parent selection. We discuss the 
properties of our compact emission line sources in a later 
paper in this series. 

The DEIMOS observations were undertaken over a se- 
ries of runs from March 2004 through April 2008. Slit 
masks were designed with position angles (PAs) within 
±30° of the measured PA in order to minimize tilt angle 
corrections in the reduction process. We used the 1200 1 
mm -1 grating blazed at 7500 A with l"slits (with excep- 
tion of 7 galaxies observed with the 600 1 mm -1 grating 
blazed at 7500 A). In this configuration, we achieved a 
spectral resolution of 1.7 A corresponding to a velocity 
accuracy of 30 km sec -1 . All DEIMOS data were reduced 
using the automated SPEC2D pipelinqj d eveloped by the 
DEEP2 survey. The spec Id packagtfjt Davis et al.|2003 l 
was used to extract ID spectra from the rectified 2D 
spectra produced by SPEC2D. The combination of ID 
and 2D spectra were analyzed using the ZSPEC soft ware, 



Coil 



axy 



also develo ped by DEEP2 (Faber et al. m prep, 
et al. ( 2004[ )), which fits a linear combination of ga' 
emission line, and stellar template spectra to each spec 
trum and allows the user to select the best-fitting tem- 
plate, thus determining the spectroscopic redshift. This 
redshift was used as the initial guess for the systematic 
velocity in our rotation curve models discussed below. 

6 http://astro.berkeley.edu/ cooper/deep/spec2d/ 

7 Based on http://spectro.princeton.edu/idlspec2d_install.html 



Because we are interested in the 2D segments of spe- 
cific emission lines in the spectra, particular care was 
taken to separate the reductions of a given mask ob- 
served at different hour angles (and therefore different 
parallactic angles), taken on different runs, or observed 
under varying conditions. A scheme (discussed below) 
was developed to coadd 2D spectra of the same target 
from multiple observation sets prior to further analysis. 
Because the spatial position of an emission line can vary 
from one observation to the next as a function of wave- 
length (by ~(f.'2), we choose to extract «100 A wave- 
length "cutouts" around key emission lines of interest 
(the [Oil] 3727 A doublet, H/3 the [OIII] 4959, 5007 A 
doublet and Ha) for our rotation curve study, based on 
the redshift determined by the ZSPEC analysis. 

Given a set of cutouts for the observation sets of 
a particular galaxy and emission line, we constructed 
optimally-weighted coadditions, with weights based on 
the signal-to-noise (S/N) and seeing full-width-half-max 
(FWHM) measured from alignment stars on the corre- 
sponding slit masks. Typically 5-6 alignment stars were 
included on each mask. Relatively sky-free windows 
(with AA < 500 A) were selected on both the blue and 
red sides of each alignment star spectrum. The stellar 
flux in these windows was weighted by the inverse vari- 
ance (as output from the reduction pipeline) and col- 
lapsed along the wavelength direction to obtain a stel- 
lar profile fitted with a Gaussian. The width and peak 
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TABLE 1 

Disk sample: Note- See Fig[T]for histogram and best 
gaussian fits of stellar mass and scale radii distribution 



Line profile 
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129 


10.11±0.05 


0.60 


2.68±0.09 


1.02 


Compact 


59 


10.21±0.10 


0.73 


2.09±0.13 


1.01 


Passive 


48 


10.69±0.09 


0.62 


2.24±0.16 


1.08 


TOTAL 


236 



















a best-fit gaussian centroid of log stellar mass 
h FWHM of best-fit Gaussian of log stellar mass 
c best-fit gaussian centroid of scale radii in kpc 
d FWHM of best-fit Gaussian of scale radii in kpc 
e possible to fit rotation curves 

were used to estimate the S/N and FWHM for the blue 
and red components of each alignment star. The average 
FWHM across a mask and the average star-by-star ratio 
of S /N values provide the seeing FWHM and relative S /N 
for that observation set. The typical spread in FWHM 
among stars on a given mask is 0.03 arcsec. The weight of 
each observation set was then given by w — s/ f 2 , where 
s is S/N and / is the FWHM, and these were normalised 
by the weights of all the coadded masks. 

All emission line cutouts were inspected by eye and 
occasionally rejected if the region extended beyond the 
detector area or if there was an artifact in the data that 
could interfere with the line of interest. To perform the 
coaddition, each 2D cutout was first rectified to a reg- 
ular grid in wavelength and spatial position using the 
2D wavelength solution output by the SPEC2D pipeline. 
We located the peak of the continuum along the spatial 
axis by collapsing the 2D cutout in the wavelength di- 
rection, initially masking out the emission feature. The 
collapsed profile was fit by a Gaussian with the result- 
ing peak taken as the galaxy center. In rare cases, the 
continuum was so weak that a position could not be de- 
termined without including flux from the emission line 
itself. The final centering of each cutout was verified 
by eye, and the cutouts were coadded with appropriate 
weighting after alignment in both wavelength and central 
continuum position. 

The seeing varied from / .'6 to 1'.' 2 over the various 
observing runs, so whenever seeing measurements are 
needed in the dynamical analysis (§31), we adopt the aver- 
age value measured from the alignment stars in the final 
coadded observations. This measurement is preferred to 
that based on a photometric image, since it not only 
refers to data integrated over the entire exposure time of 
the spectra, but we also account for systematics in the 
observing and co-addition process (which use different 
weights for different exposures). 



2.3. Photometric data and bulge-disk decomposition 

By selecting our sample within the GOODS North and 
South fields, we ensure excellent quality multi-color data 
for all our galaxies from the Hubble Space Telescope 
Advanced Camera for Surveys (HST ACS) ( jGiavalisco] 
et al. 2004). This provides valuable structural informa- 
tion that can be used for translating the observed ro- 
tation curves into physically-based properties. To the 
imaging in four bands from HST (-B435, Veo6, 2775, and 
Z850 -bands), we add ground-based K-band data in order 
to se cure stellar mass estimates based on SED fitting (see 

pa. 



A key requirement for our analysis is the inclination, 
PA and effective radius of each galaxy. We also need to 
separate, where possible, the disk light from bulge con- 
tamination. We derive these quantities from the HST 
imaging using the galfit3 ( |Peng|[20T0l ) least squares 
elliptical-fitting method. For each galaxy we extracted 
a 9V03 x 9'.'03 (301 x 301 pixels) postage stamp centered 
on the object. Neighbors were individually masked out 
to eliminate confusion. We first fit an exponential disk 
component plus a dc Vaucouleurs' bulge profile to every 
galaxy. Those galaxies which yielded unphysical solu- 
tions were re-fit with a single Sersic profile component, 
where the Sersic index (n) was allowed to vary. Such 
cases generally represent disk galaxies which are bulge- 
less and / or more clumpy and irregular than regular well- 
formed spirals. ~60% of our galaxies were fit using a 
1-component n- varying Sersic profile fit, and ~40% were 
adequately fit with a two-component bulge and disk so- 
lution. This mixture gives some indication of the mor- 
phological distribution of our sam ple in dicating that less 
than half are well-formed spirals ({ 2.1 ). Disk sizes, incli- 
nations and PAs were taken from best-fit disk component 
if more than one component can be fit. To track possi- 
ble biases in the disk-bulge decomposition we will later 
flag those galaxies for which significant bulge components 
were measured. 

We ran GALFIT using HST data in all four bands 
(Bp435W, Vfgoqw, IF775W, an d zps50Lp)- The disk 
scale radii are consistent among the bands indicating no 
significant redshift-dependent bias (or morphological k- 
correction) within the sample. In order to maximize the 
signal-to-noise, hereafter we used the GALFIT results from 
the zps50LP band. To achieve convergence on the gal- 
fit parameters and to assess any systematic uncertainty 
in the fitting technique, we ran a Monte Carlo analy- 
sis (N=1000) where we varied the initial guess of each 
parameter (magnitude, effective radius, b/a for inclina- 
tion, position angle, and sky) from one adopted by th e 
GOODS S Extractor results ( |Giavalisco et al.| ( |2004[ )). 
We found the parameter output distributions were much 
narrower than the input distributions thereby demon- 
strating convergence. Final parameter uncertainties from 
the Monte Carlo distributions are better than 5% on av- 
erage, and we add these uncertainties in quadrature to 
the observational error. 

2.4. Stellar masses 

Reliable stellar masses are an essential component of 
constructing a baryonic Tully-Fisher relation. We take 
o ur stellar mass estim ates from work initially presented 
in |Bundy et al. J 2005 ) , followed by the analysis presented 
in jBundy et al.| ( |2009 |. Further details can be found in 
those papers. 

Briefly, stellar masses are derived using a matched cat- 
alog of multi-band ACS and ground-based K 5 photome- 
try. The essential near-infrared data was taken with the 
MOIRCS imager on the Subaru telescope for GOODS-N 
(Bundy et al. 2009[ ) and the ISAAC instrum ent on the 
ESO VLT for GOODS-S ( |Retzlaff et al |2010[ ). The final 
matched catalog is substantially complete to a limiting 
magnitude of i^AB=23.8, deeper than our spectroscopic 
limit. 

A Bayesian code fits the spectral energy distribution 
(SED) derived from 2 arcsec diameter ACS and K s pho- 
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tometry adopting the best spe ctroscopic redshift and this 
SED is compared to a grid of Bruzual & Chariot ( 2003 ) 
models that span a range of metallicities, star formation 
histories, ages and dust content. The stellar mass is es- 
timated by multiplying the derived K-band mass/light 
ratio by the observed K-band luminosity derived from 
the MAG_AUT0 total K ron magnitude d etermined by SEx- 
tractor. We assume a Chabrier ( 2003 ) initial mass func- 
tion (IMF). The probability for each fit is marginalized 
over the grid of models giving a stellar mass posterior 
distribution function, the median of which is the cata- 
logued value. At the magnitudes probed in this survey, 
the uncertainties inferred from the width of these poste- 
rior functions is less than 0.1 dex. Including systematic 



errors (see Bundy et al. ( 2005 ) for a full discussion) , we 
determine that the stellar masses are reliable to better 
than 0.2 dex, modulo uncertainties arising from the IMF 
normalization. 

In order to construct a self-consistent TF relation, we 
apply an aperture correction to the total stellar mass esti- 
mates (i.e., for a given fiducial radius, an aperture stellar 
mass at that radius compared to the velocity measured 
at that fiducial radius). We extract equivalent Kron ra- 
dius aperture fluxes on the zp850LP band data to get the 
flux equivalent to that used in the total stellar mass esti- 
mates. We then take a scaled-down aperture flux within 
the given fiducial radius, and compare this to the Kron 
radii flux thus deriving a total-to-enclosed flux ratio. As- 
suming the ZF850LP-band and iC-band are roughly equiv- 
alent stellar tracers, we can use this ratio to estimate the 
enclosed stellar mass. This approach maximizes the util- 
ity of the HST ZF&50LP images which have much better 
resolution than our ground-based K-band data, thereby 
giving us the resolved mass distribution throughout the 
disk to match the detail seen in our rotation curves. 

3. DYNAMICAL ANALYSIS 

We now turn to the questions of how to extract re- 
liable rotation curves from our 2-D spectroscopic data 
and how to interpret those curves in terms of a fiducial 
velocity measurement that can be used in the various 
Tully-Fisher relations we will example. To fully exploit 
our extended integrations, we re-examine the rotation 
curve model and consider carefully how to define a self- 
consistent fiducial velocity that can be robustly measured 
within our data. We quantify improvements in our data 
by comparing our velocities (extracted from 6-8 hours of 
integrated exposure time) with those deduced from ~1 



studies flVogt et al.||1997[ 


Conselice et al. 


2005 Weiner 


et al.|2006b Kassin et a,l.\'z 


007 


) . We also examine the re- 



the first studies with IFU spectrographs. In what fol- 
lows, our analysis is based on t he 12 9 galaxies for which 
extended emission is observed (|2.2[). 



3.1. Rotation curve model 

Optimally fitting rotation curves presents a variety 
of challenges that become more difficult as redshift in- 
creases. Foremost, we seek a functional form which rep- 
resents the bulk of the observed emission line shapes and 
has some physical basis. Secondly, we must aim to char- 
acterize this functional form with a fiducial velocity that 
is reliably detected across the sample, preferably without 
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Fig. 2. — A functional arctan rotation curve with various charac- 
teristic radii (discussed in the text) shown as dimensionless factors 
of the disk scale length (r^ 1 ). The example shown has a turnover 
radius, rt ~ 0.4r s , and an asymptotic velocity, V a ~ 130 kms -1 , 
typical of many galaxies in our sample. The region where the rota- 
tion curve turns over is extensive and the so-called turnover radius, 
rt , does not necessarily indicate the most appropriate center of this 
region for the arctan function. The histogram below shows the ex- 
tent to which we can reliably trace emission lines in the spectra 
of our sample. The frequently used V m ax (maximum measured 
velocity) is not equivalent to the asymptotic value, V a , a math- 
ematical extrapolation and typically not reached in the observed 
rotation curve. Estimates of V a can depend critically on how well- 
constrained are the central region and turn over of the rotation 
curve. As ~90% of our galaxies are traced to 2.2r s (and all to at 
least 1.0r 3 ), in our TF relations we opt to use V2.2, the velocity at 
r which minimizes uncertainties arising from extrapolation (see 

extrapolation to radii where there is no data. Given our 
extended integrations, we have considered carefully the 
optimum selection of this characteristic velocity. The 
challenges can be appreciated by considering Figure [2] 
where we show various characteristic velocities in the 
contex t of the frequent ly- used arctan model of a rotation 
curve ( Courtea"u||1997 l as well as the extent to which we 
can trace emission lines for our sample. 

Several functional forms have been discussed in the 
literature, such as the "multi-parameter function" in 



Fersic et a 



Courteaul J1997 1 and the "universal rotation curve" of 
( |l996 1. The simplest model flexible enough 



to fit most rotation curves is the empirically-motivated 
arctan function (see Fig. [2]), which we adopt here, viz: 



2 , r — TYi . 

V = V + -V a arctan 

7r r t 



(1) 



where Vb is the central or systematic velocity, ro is the 
dynamic centre, V a is the asymptotic velocity, and r t 
is the turnover radius, which is a transitional point be- 



Courteau 


1997 


Willick 


1999 



bulge-dominated rotation curves around the turnover ra- 
dius, but we typically do not observe this feature in our 
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Fig. 3. — An example of our fitting procedure for a galaxy with asymmetric emission. The left panel shows the observed 2-D spectral data 
in the region of Ha, and the right panel is the equivalent section of the constructed model spectrum. Both show a grey scale representing 
flux as a function of velocity with respect to position. Orange circles show the trace centroids of the emission line (bars along the positional 
dimension indicate bin size, not error), and the green triangles show the equivalent centroids for the constructed model. The T2.2 radius is 
marked by vertical dotted lines and the extracted V2.2 velocity is marked by horizontal solid lines. These lines of fiducial measurements do 
not cross directly through the t race because they are extracted from the best-fit model's intrinsic arctan function, prior to distortion by 
seeing and dispersion. See §3.2| for more details. 

rotation curves flatten is visually confirmed in most of 
our galaxies. TF velocities based upon V2.2 result in the 
smallest internal scatter and provide the best match to 
radio (21 cm) line widths for local galaxies ( jCourteau 
1997[). This fiduc ial velocity has been adopted by|Uut 



sample. 

Past studies employed as TF velocities the maximum 
measured velocity along the di sk, V m ax, or the asymp- 

^ T996j [19971 
e disadvan- 



totic velocity from Eq. Ill V a flVogt et al. lITi 
Weiner et al.|[2006b Kassin et al. 20071). The 



Weiner et aL]|2006b| [Kassin et al.||2007| ). T 

tage with V max is that it is not measured at a consistent 
location in the variety of disks observed. In terms of the 
disk scale radius, we can see a range of a factor ~5 or so 
in the associated radius. 
Some studies (i.e 



2006 Weiner et al. 



Flores et al. 

|2006a[ ) advocate the use ot V C irc, the circular radial ve 
locity, driven by either the disk or the halo, but in 
nearly all cases this is assumed to be the V a of the 
widely adopted arctan model. While the arctan function 
closely matches the observed extent of rotation curves, 
the mathematically-extrapolated asymptote is unwar- 
ranted as it is not observed in typical datasets, even for 
our extended integrations. Small offsets in the extrap- 
olated velocity curve will lead to large changes towards 
the asymptotic limit, especially when emission is not de- 
tected past the flattened part of the rotation curve, and 
emission line dispersion and seeing distort the terminal 
emission. Under these circumstances, the majority of the 
l V C i r c velocities remain as extrapolations. 

We thus seek a physically motivated fiducial radius to 
which emission can be detected across the variety of disks 
seen over our wide redshift range. We have adopted the 
modeled velocity at 2.2 times the exponential disk scale 
length, r 22 , which we will call ^2.2- This has a good 
physical basis as it corresponds to the location of peak 



rotational amplitude for a pure exponential disk ( Free 

T9871 ICourteau & Rix 



man 



[19701 |Binney fc Tremaine| [l 
. Although few disks are lik 



ton et al. ( |2010a[ ) in determining the kinematic connec- 
tion between late-type galaxies and dark matter halos. 
As shown in Figure [2] we trace the velocity field with at 
least one optical emission line to r 2 . 2 for ^90% of our 
sample (Fig. [2|. 

3.2. Rotation curve fitting procedure 

We now discuss how we fit the emission line data to 
obtain an accurate velocity measurement at r 2 . 2 for each 
galaxy disk using the arctan model. In the next section 
we discuss various inclination corrections that must be 
made following the fitting procedure. A simple outline 
of our fitting code is as follows: 

• Trace the observed 2D emission line in wavelength 
as a function of spatial position, 

• Construct a model 2D spectrum with an arctan- 
shapcd emission line profile, implementing features 
measured from step 1 (e.g., position-dependent 
emission brightness profile, dispersion, and seeing), 

• Trace the model spectrum, varying the arctan 
model parameters V a and r t (as well as the dynam- 
ical center and a seeing-dispersion beam correction 
factor) until the model trace optimally fits that of 
the data, 



1997 ) . Although few disks are likely to be perfectly ex- 
ponential, the r 2 . 2 approximation as the point at which 



• Compute V2.2 at r 2 . 2 
els. 



from the best-fit arctan mod- 
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Fig. 4. — Examples of data from the survey. For each source we show an HST ACS composite (B, V, i,z) photometry postage stamp 
rotated for convenience so the overlaid 1 arcsec slit level is horizontal, the Keck DEIMOS 2-D emission line cutout, the corresponding best 
model and residual of model minus data. The extent of T2-2, where V2.2 is measured is also overlaid. The triptych second to the bottom 
shows the same emission line as the one above but with only 1 hour spectroscopic integration. 



In detail, our procedure is as follows: 

Step 1: We begin by subtracting the galaxy contin- 
uum on the rectified 2D spectral frame, leaving only the 
emission-line profile. This is performed by linearly inter- 
polating in wavelength across the emission line region, 
treating each spatial row independently. 

Step 2: We next trace each emission line as a function 
of spatial position. This procedure returns the central 
wavelength, line width, and peak flux in each spatial bin. 
We fit two half-Gaussians to the line profile in the wave- 
length direction, and use an adaptive binning procedure 
in the spatial direction to ensure a consistent signal-to- 
noise (S/N). The blurring effect of seeing in the spatial 
direction, and of dispersion in the wavelength direction, 
can affect the emission line in a particular spatial bin 
by mixing flux from neighboring bins. We use two half- 



Gaussians with differing sigmas joined at the same peak 
to account for these effects, ensuring that the position of 
the peak flux is always traced. 

The trace terminates when the emission is no longer 
detected above the local noise level. All traces are visu- 
ally inspected to ensure that spurious reduction artefacts 
do not interfere with the fit. The procedure produces an 
array of wavelength values as a function of spatial posi- 
tion (see Fig. |3| , as well as the emission peak flux. 

We also experimented with different forms for the emis- 
sion line brightness profile along the slit, for example the 
collapsed light through a synthetic slit estimated from 
the HST image, or an exponential profile. However, 
traces using these profiles resulted in larger residuals. 

Step 3: We next construct a model 2D spectrum from 
the rotation curve of eqn. 0, which requires V a and 
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Fig. 5. — As Figure[4]but for more disturbed, visually asymmetric disks, some of which appear to be undergoing minor mergers. Although 
the morphologies are more irregular, we usually succeed in fitting an appropriate arctan-based model. 



parameters for the arctan functional form. We guess ini- 
tial values to approximately match the flat portion of 
the outer trace of the observed data. We lay this model 
rotation curve on a 2D grid, adjusting the normalisa- 
tion (brightness) and line width in each spatial bin to 
match that of the data (i.e., using the parameters from 
the two half-Gaussians in each bin) . We finally conv olve 



the resulting rotation curve with the local seeing (see [2.2 



for details), and the resulting peak model brightness is 
re-normalized, bin by bin, to that observed. Ideally we 
would start with the unblurred emission-brightness spa- 
tial profile, but there is a degeneracy in blurring by see- 
ing and dispersion, which is particularly troublesome for 
galaxies with irregular emission line brightness profiles. 

Because we have measured the position-dependent dis- 
persion from the trace (which implicitly contains the ef- 
fect of the seeing) , as well as a separate measurement of 
the seeing from the alignment stars, we can attempt to 



break this degeneracy by fitting for a multiplicative fac- 
tor. This correction factor is used to multiply the disper- 
sion implemented in the model, and it is always less than 
or equal to 1. This allows us to approach the intrinsic, de- 
convolved rotation curve. This method is more success- 
ful than alternative approaches we attempted using both 
synthetic and observed data, for example adopting the 
collapsed light profile from the most appropriate broad 
band HST image as observed through the slit, or simple 
functional forms of emission brightness profiles (constant, 
linear, exponential). Indeed, we find many emission lines 
that are brighter at the disk edge than towards the cen- 
ter, as well as asymmetric emission distributions which 
do not match the broad-band flux distributions that en- 
ter the slit. 

Step 4 : To select the best model, we use a robust 
non- linear least squar es fitting algorithm based on MP- 
fit ( Markwardt|2009 1 . The arctan input parameters, V a 
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and Tt are varied as well as the position of the dynamical 
center and the seeing-dispersion beam correction factor. 
Chi-squares are calculated between the model trace and 
the data trace, rather than on the entire 2D frame. Not 
only is this method much faster, but we find it results in 
smaller residuals, simply because we are focusing the fit 
on the information that is most important to the shape 
of the rotation curve. We propagate the error from the 
measured input parameters (emission profile, seeing, and 
dispersion) by using a Monte Carlo approach, simultane- 
ously altering input parameters from random Gaussian 
distributions with widths matching the error in the ob- 
served parameters for each galaxy. After 100 iterations, 
we add the uncertainty from the resulting distribution of 
output velocities to the formal fitting errors. 

The velocity at r 2 .2 can then be calculated using the 
best-fit arctan model and the r adius measured from the 
ZF850LP band HST data (see (2.3). All the observed 
emission lines in a given spectrum are treated indepen- 
dently in the above steps. These independent measure- 
ments of V2.2 are found to be consistent in > 90% of our 
galaxies. We combine the emission line fits for the same 
galaxy into a final weighted average of V2.2 and include 
the error in r 2 .2 in the final uncertainty. 

Examples of our spectral data and respective best- 
fit model rotation curves can be found in Fig. [4] and 
Fig. [5] Note that these models represent the minimum 
chi-square best-fit on the 1-D trace of the 2-D models 
as described above, not a direct 2-D chi-square best- 
fit. Fig. [4] shows six examples of disks which appear 
mostly morphologically-regular in their HST multi-filter 
photometry, whereas Fig. p] displays six galaxies which 
are likely more disturbed or irregular based on their HST 
images. In both figure sets, we observe fairly regular 
arctan-shaped kinematics, although places of high dis- 
persion and brightness in the emission line tend to coin- 
cide with regions of the galaxy in the HST image that 
appear to be internally or externally disturbed. For the 
5th galaxy in both figures, we show the best-fit model 
rotation curve after only 1 hour of integration time for a 
comparison to the full 6-8 hours of integration shown in 
the panel just above. 

3.3. Inclination and PA offset corrections 

We now correct our V2.2 measurements for the effects 
of disk inclination and any misalignment between the 
position angle (PA) of the DEIMOS slit and the major 
axis of the galaxy as determined from GALFIT. 

Adopting the convention i = 0° for face-on and i = 90° 
for edge-on disks, the inclination correction is: 



v , 



(sini) ' 



= cos 



! (b/a) 



1-98 



(2) 
(3) 



where go— 0.1-0.2 represent s the intrinsic flattening ra- 
tio of an edge-on galaxy dHaynes k, Giovanelli 1984 



Courteau|[l996) |Tully et al'||1998| ). Although the pre 



cise value depends on morphology, the uncertainty leads 
to changes in t he final velocity measurement on the or- 
der of 1 kms" 1 (Pizagno et al. 2005| |Haynes fc Giovanelli 



19841 



Pizagno et al 
For the PA offset 



=0.19 for all systems, similar to 



We assumed qo= 
p005]. 

we determine a correction: 



Vol, 



cos (APA) ' 



(4) 



from simple trigonometry accounting only for the mis- 
aligned slit component of the true major axis. Only 10 
galaxies in our final sample have velocity corrections aris- 
ing from a PA offset greater than 10%, and no PA offset 
exceeds 45°. 

We apply these two corrections to produce a catalog 
of corrected V2.2 measurements, where the effects of see- 
ing, velocity dispersion and the emission line brightness 
profi le are accounted for in the model. As explained in 
(2.4 to compare like with like in our TF relations, we 



match V2.2 to an estimate of the stellar mass within r2. 2 - 
In order to compare this enclosed relation to a more fa- 
miliar construction of the TF relation used in previous 
studies, we also consider a total relation, which com- 
pares the stellar mass estimated from the Kron radius 
with the velocity associated with the optical radius rjj, 
defined for an exponential profile as one enclosing ^83% 
of the light, or 3.2r s . This measurement of V3.2 or Vr> we 
find to be comparable to the V C i rc or V max of previous 
intermediate-redshift studies, however our choice of rjj 



is more consistent between disks than the 
those studies. 



used in 



3.4. Demonstrating the advantages of extended 
integrations 

A major advance in our survey is the use of extended 
integration times (typically 6-8 hours) on sources with 
apparent magnitudes (zp850LP < 22.5) that have typ- 
ically been observed for ~1 hour exposure time ( Vogt 
et al.||1997| IConselice et al.| |2005| |Weiner et aLp OOeb; 
Kassin et al.| 2007 ). In addition to ensuring our emission 
fines are traced to r 2 .2 (Figure[2|, this leads to improved 
signal/noise at all radii. We can demonstrate the effect 
this has on our derived velocities by using our model fit- 
ting code on a subset of our data taken with a 1 hour 
integration time. Comparing the results with the 6 hour 
integrations in Figure [6j we can draw two important con- 
clusions. 

Firstly, while most V2.2 measurements are consistent 
given the error bars, the 1 hour data has significantly 
larger error bars and there is a systematic offset whereby 
the 1 hour analysis produces on average a 13% higher 
measurement of V2.2 than the 6 hour data. This suggests 
that if errors are not properly accounted for, studies us- 
ing a similar modeling procedure but with less integra- 
tion time may falsely detect evolutionary signals as a 
result of over-estimated velocities due to the decline in 
signal-to-noise. 

Secondly, as expected, the 6 hour data enables us to 
probe largely beyond T2.21 whereas this is not the case for 
the 1 hour subset. The right hand panel in Fig. [6] shows 
the gains made in detected emission extent in factors of 
scale radius. We trace, on average, 30% further along 
the disk with the 6 hour data than the 1 hour data. In 
fact, modeling an arctan function fails for one third of 
the 1 hour sample due to the low signal-to-noise, and 
missing segments of the emission create a false deviation 
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Fig. 6. — A comparison of the modeled V2.2 values and maximum measured emission line extent as derived from our total exposure time 
(6 hours in these examples) to those derived from a 1 hour exposure extracted as a subset of our data. Left: the agreement between V2.2 
values is reasonable given the errors, however there is a systematic offset and the errors are significantly larger for the 1 hour subset. Right: 
almost every galaxy reveals significantly more extended emission in the 6 hr data. The shading marks the r < r^.i region; if the emission 
does not extend beyond this region, V2.2 must be extrapolated, as is the case for several lhr galaxies. Note also that a third of the rotation 
curves cannot be adequately traced with only 1 hour of integration (circled points). 



from the arctan shape. We are confident that we have 
successfully converged on necessary integration time in 
our 6 hour data because emission is traced beyond r-2.2 
in ^90% of our disks, ensuring we can accurately measure 

V2.2- 

3.5. Comparison to previous work 

Our rotation curve fitting method was constructed to 
work on data with extended exposure times, so it is in- 
teresting to compare to velocities for the same galax- 
ies derived using alternative techniques. The greatest 
overl ap can be found using the TKRS/G OODS sam- 



V 300 



2004 Giavalisco et al. 2004). Fig- 



shows 35 of our galaxies in common with in the 



pie (Wirth et al. 
ure [7] 

TKRS/GOODS sample for which there are 2D spectral 
fits for V ro t (equivalent to V a and does not include a 
correction for inclination), and all of our sample over- 
laps with their ID linewidth measurements ( Weiner et al. 
2006a|b[ ). We compare the velocities we measure for these 
galaxies without the inclination correction applied for a 
most direct comparison. 

On average, the mean V ro t for the TKRS 2D data is 
68% of the mean V2.2 for our study. For the Id linewidth 
data, the TKRS measurements have a mean that is 75% 
of the V2.2 mean in our study. Our current data have sig- 
nificantly smaller error-bars as would be expected given 
the longer integration times: the median fractional error 
bar is 8.6% for this study, and 59.1% for T KRS. 

The three galaxies in com mon with the |Flores et al.| 
(2006); Puech et al. (2008) studies have consistent ve- 
locity measurements within the error bars when compar- 
ing only the equivalent slit area from our study to their 
modeled IFU data. One galaxy has an inconsistent TF 
velocity measurement since a higher V max is found on 
either side of the slit area on the full modeled IFU ve- 
locity field. Little can be concluded from a comparison 
of three objects, however in terms of sample selection, it 
is encouraging that the three shared objects between our 
studies all belong to the sub-class of Complex Kinemat- 
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Fig. 7. — A comparison of the velocities for 35 galaxies tha t over- 
lap be tween thi s study and the TKRS/GOODS study of |Wirth| 
let al.|(|2004k & iGiavalisco et al.l d2004l. In black are shown the 
The 'l'kRS/GOUUS '2D spectral fits lor V ro t (equivalent to V a ), 
and in gray are the TKRS/GOODS ID linewidth measurements. 
Inclination corrections are not applied to any of the velocity mea- 
surements plotted here. 



in the Flores et al. (2006) kinematical classification 



scheme, consisting of objects with the most irregular, pe 
culiar observed kinematics. 

4. RESULTS 

We now reach the primary aim of our paper: to present 
the redshift-dependent stellar mass Tully-Fisher (M*- 
TF) relation over 0.2 < z < 1.3. We will explore 
the quantitative improvement, which we have realized 
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z[0.2-0.5] z[0.5-0.6] z[0.6-1.25] 
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Fig. 8. — Redshift-dependent stellar mass Tully-Fis her (M^-TF) relations usin g the enclosed (top pa nels) and total (bottom panels) mass 
estimates. Dashed lines refer to the local relations oflBell & de Jongl (12001 ft andlPizagno ct al. (2005). Dot-dashed lines refer to the fit to 
the lowest redshift bin in the current dataset. To fit fixed slopes between redshift bins, we tit relations using a Monte Carlo distribution of 
slopes from the best-fit free slope of the entire sample. The resulting mean slope is shown as a solid red line with red dotted lines showing 
the la error in the zero-point (in M*/Mq dex). Using the velocity and enclosed mass at the ri.i aperture reduces both the intrinsic scatter 
((Tint) and the rms of the relationship in each redshift bin. 



TABLE 2 

Stellar Mass Tully-Fisher Relations 



z range 




N 


r\ 


¥\ 












median 






Enclosed: M,(r 2 . 2 ) v» V(r- 22 ): 


0.2<z<1.3 


0.64 


129 


1.718±0.415 


3.869±0.193 


0.058 


0.022 




0.083 


0.224 


0.091 




0.323 


0.2<z<0.5 


0.37 


43 


1.755±0.043 


" fixed 


0.055 


0.035 




0.077 


0.211 


0.081 




0.297 


0.5<z<0.8 


0.62 


49 


1.684±0.049 




0.067 


0.045 




0.086 


0.257 


0.118 




0.334 


0.8<z<1.3 


0.96 


37 


1.720±0.053 




0.052 


0.062 




0.086 


0.202 


0.098 




0.331 


Total: M* vs V(r D ) 


0.2<z<1.3 


0.64 


129 


1.926±0.476 


3.783±0.223 


0.070 


0.023 




0.093 


0.266 


0.087 




0.353 


0.2<z<0.5 


0.37 


43 


1.969±0.050 




0.069 


0.037 




0.088 


0.262 


0.083 




0.332 


0.5<z<0.8 


0.62 


49 


1.902±0.050 


" fixed 


0.072 


0.057 




0.093 


0.272 


0.110 




0.351 


0.8<z<1.3 


0.96 


37 


1.900±0.063 




0.073 


0.064 




0.099 


0.277 


0.084 




0.375 



a best-fit y-intercept in M*/Mq dex assuming scatter in V/kms -1 dex 

b slope assuming scatter in V/kms -1 dex 

c internal scatter in V/kms -1 dex 

d median velocity error in V/kms -1 dex 

c total scatter in V/kms -1 dex 

'internal scatter in M,/Mq dex 

g median stellar mass error in M*/Mq dex 

h total scatter in M»/Mq dex 
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through our extended integrations and subsequent ro- 
tation curve modeling, with respect to earlier work in 
a number of ways. Foremost, the scatter around the 
redshift-dependent relations will provide a good indica- 
tion of our progress. We will examine relations using 
masses derived within our chosen fiducial 2.2 scale radius, 
J*2.2j the so-called enclosed relations alo ngsi de those for 
the mo re traditional total relations (see {2.4 and the end 



of £3.1 for more details). It is also convenient to exam- 
ine and discuss the i?-band luminosity-based TF relation 
from our survey as the literature contains many estimates 
of this scaling relation and previous studies have claimed 
evolution, de spite large scatter and possible incomplete- 



ness biases (£4.2 ) 



4.1. The redshift-dependent stellar mass TF relation 

We begin by considering the case for evolution in the 
stellar mass (M*)-TF relation. The results are illustrated 
in Figure [8] and listed in Table [2] Redshift bins were 
selected to ensure nearly uniform samples over our total 
redshift range. Shifting the boundaries of these bins by 
modest amounts does not change the overall conclusions 
we present below. 

To fit a linear regression to our data we adopt a least- 
squares approach which incorporates a measurement of 
the intrinsic scatter o"; n t, which is added in quadrature 
to the velocity dimension. We fit a zero-pointed line: 



log (M,) = [a+ 6 log (V2.2)] - log (M ), 



(5) 



where Mq = 10 M©, and while we plot the relation in 
the familiar way with velocity on the x-axis and a y- 
intercept given in terms of stellar mass, we treat velocity 
as the dependent variable in the linear regression. Fit- 
ting linear regressions with stellar mass as the dependent 
variable leads to fits which s uffer much more from the ef- 
fects of incompleteness b ias ( |Bamford et al.|2006 Weiner 
eTaLl|2006a| [KelTy|[2U07| . 

To fit the data to Eq. [5j we adapted a code which 
takes into account errors in both the ordinate and ab- 
scissa. We first fit an unrestricted slope to the entire, 
un-binned log V - M* data set, and then fit the individ- 
ual redshift-binned relations with a distribution of slopes 
drawn from a Monte Carlo (N=100) Gaussian distribu- 
tion of slopes centered on the slope found for the full 
sample. We undertake this exercise for both the enclosed 
mass (i.e. that associated with our fiducial radius, ^.2), 
and the total mass. We also experiment with simpler 
approaches to linear regressions, which result in similar 
trends but naively smaller errors, so we present here our 
most robust results with the broadest consideration of 
the uncertainties. 

For a local comparison, we consider the M*-T F relation 
derived from i^-band luminosities publishe d by|Bell fc de 



Jon, 



pOOTT) and the M*-TF relation fromlPizagno et al. 



2005 ) , the latter of which is based on velocities from ri.i 



as presented in the enclosed relation of this study. 

Before discussing possible evolution, we consider the 
derived scatter around the relations since this is a valu- 
able indication of our gains in precision. Satisfactorily, 
we find intrinsic scatters in the fits of 0.055, 0.067, 0.055 
(dex of V km s^ 1 ) at (z)=0.37, 0.62, 0.96, which are com- 
parable to tha t seen in local TF relations (i.e., ^0.049 in 
Pizagno et al. ( |2005[ )). Our analysis therefore represents 



a significant improvement on the scatter seen in earlier 
intermediate-redshift studies, for examp le an improve- 
ment of a factor of 2-3 over the study of IConselice et al.] 
(2005). We also achieve a scatter less t han t hat found 
in the M*-5o.s relation of Kassin et al. (20071) (between 
0.08 and 0.11 dex), despite the tact that they introduce 
an additional dispersion term, So. 5, which significantly 
tightens the relation from the log V - M* relation. 

Table [2] shows little room for evolution in the rela- 
tion. Fitting a straight line through the zero-points be- 
tween redshift bins, we detect a modest but statistically 
insignificant trend for a larger stellar mass at fixed ve- 
locity at lower redshift: AM*~0.037 ± 0.065 dex (la) 
from (z)~1.0 to 0.3. Fitting a stra ight line through zero- 
points between redshift bins in the Conselice et al. (2005 
study, we see a AM*~0.07± 0.19 dex (lcr) from (zj^l 
to 0.3, consistent with our result, yet more uncertain. In 
terestingly, our 1-sigma limit for evolution is consistent 



with the modest AJog M*/M© predict ed by Portinari fc 



|Sommer-Larsen| ( |2007[ ) (AM* ~ 0.1), |Somerville eTal 



( |2008p (10%~Hecreas e in V at fixed Mj w ith time), and 
only consistent with Dutton et al. (2011a) (AM* ~ 0.2) 
over a similar redshift range at the 2-sigma level. Impor- 
tantly, this result is robust to the inclusion of the local 
data in both the enclosed and total relations. We discuss 
the implications of this agreement further in £|5j 

We note that the scatter in the total stellar mass TF 
relations is increased somewhat compared to that in our 
preferred enclosed relations, most likely due to the effect 
of extrapolated velocity measurements. We trace spec- 
troscopic emission beyond rr> on ^60% of our disks as 
opposed to ^90% beyond r-i.i- 

4.2. The B-band magnitude TF relation 

W e undertake a similar analysis to that described in 
j |4.1| for the absolute B magnitude (Mb) TF relation, 
where Mb me asur ements come from the best-fit SEDs 
(described in \2A), and are aperture corrected in the 
same way as the stellar mass estimates. The TF rela- 
tions are shown in Fig. M for both total and enclosed 
luminosities and the results are listed in Table [3] Any 
difference in the redshift-dependent trends compared to 
that for stellar mass relation would indicate changes in 
the star formation rate per unit stellar mass. Some work- 
ers have claimed to see evolution in the Mg-TF relation 



( Weiner et al.|206"6a Fernandez Lorenzo et al.|2010 | and 
we aim to verify or otherwise these trends with our im- 
proved dataset . As b efore, we use the local re lations of 
Pierce fc TuhVl ( p92| and the (z)~0.3 study of |Bamford 
et al. 1 2T)06|) as comparison datasets. 



Once again, the intrinsic scatter around our redshift- 
dependent M B -TF relations, 0.424, 0.641, 0.670 in mag, 
at (z)=0.37, 0.62, 0.96, are com parable to those seen in 
the local relat ions (~0.4 mag i n [Pierce &: Tully (1992), 
and 0.3-0.5 in Verheijen (2001)), and we note a consid- 



erable improvement over pre vious intermediate-re dshift 
studies. The (z)~0.3 s tudy of i Bamford et al. (2006) and 
the (z)~0.85 study of (Jhiu et al.| ( |2007 ) have scatters 
twice a s large (~0.9 in mag), and those of Weiner et al 



(|2006bj) to z~1.2 and Fernandez Lorenzo et al. "(2010) to 
z~1.4 nave scatters ~2-3 times as large (~1.5 mag and 
~1.2 mag) respectively. 

Even though the M^-TF relation is not as tight as 
our M*-TF relation, there is evidence for a stronger evo- 
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z[0.2-0.5] z[0.5-0.B] z[0.8-1.25] 




1.8 1.8 2.0 2.2 2.4 2.6 1.6 1-8 2.0 2.2 2.4 2.6 1.6 1.8 2.0 2.2 2.4 2.6 

log(V/kms"') 

Fig. 9. — - As Fig. EH but for th e absolute B-band magnitud e (Mr) Tully-Fisher r elation. For comparison purposes, we show the local 
relation of |Pierce &; T ully (1992) and the z ~ (0.3) relation of Bamford ct al. (2006) as dashed lines. Other lines are as in Fig. [8] 



TABLE 3 

Absolute B-band Magnitude Tully-Fisher Relations 
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13 















median 




~M3 


Enclosed: 


C T "2.23 vs V ( T '2.2) ; 


























0.2<z<1.3 




0.64 


129 


-3.589±0.893 


-7.546±0.578 


0.081 


0.022 




0.127 


0.612 


0.291 




0.956 


0.2<z<0.5 




0.37 


43 


-3.260±0.104 




0.056 


0.035 




0.094 


0.425 


0.245 




0.711 


0.5<z<0.8 




0.62 


49 


-3.703±0.133 




0.085 


0.045 




0.130 


0.641 


0.364 




0.979 


0.8<z<1.3 




0.96 


37 


-3.906±0.162 




0.089 


0.062 




0.134 


0.670 


0.299 




1.011 


Total: M B vs 


V(r D ) 


























0.2<z<1.3 




0.64 


129 


-3.413±1.037 


-7.754±0.657 


0.091 


0.023 




0.130 


0.706 


0.250 




1.008 


0.2<z<0.5 




0.37 


43 


-3.127±0.114 




0.065 


0.037 




0.099 


0.505 


0.258 




0.771 


0.5<z<0.8 




0.62 


49 


-3.545±0.144 




0.094 


0.057 




0.132 


0.731 


0.364 




1.027 


0.8<z<1.3 




0.96 


37 


-3.624±0.185 




0.109 


0.064 




0.144 


0.845 


0.272 




1.114 



a best-fit y-intercept in mag assuming scatter in V/kms -1 dex 

b slope assuming scatter in V/kms -1 dex 

c internal scatter in V/kms~ dex 

d median velocity error in V/kms^ 1 dex 

e total scatter in V/kms~ 1 dex 

'internal scatter in mag 

g median B-band magnitude error 

h total scatter in mag 
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lution in the M^-TF relation than in the M*-TF rela- 
tion, as expected from the well-es tablished increase in 
disk star formation rate to z ~ 1 (Bundy et al. 2005). 
Fitting a linear regression through the zero-points be- 
tween redshift bins of the enclosed Mg-TF relation we 
find AMg~0.85iO.28 mag evolution in the relation from 
(z)^1.0 to 0.3. We can check if this result is affected 
by a Malmquist bias (given the distribution of luminosi- 
ties is significantly different between redshift intervals), 
by comparing subsets with similar luminosities and stel- 
lar masses, and the t rends do not substantially change. 
Weiner et al.| ( |2006b ) find a consistent trend of ~1.0-1.5 
mag evolution from a similar redshift range of (z)~l.2 to 
0.4, but with more uncertainty. Our results a re consis- 
tent with the evolution s hown in the models of Po rtinarH 
fc Sommer-Larsen| ( |2007[ ): AM B ~ 0.85 from z ~ 1. We 
can understand the different evolutionary trends in the 
Mg-TF and M„,-TF relations by examining the redshift- 
dependent correlation between Mb and Af*. To first or- 
der, as expected, the difference is explained by the in- 
crease in the B band luminosity per unit stellar mass 
with redshift ( |Lilly et al.||1996| |Madau et al.||1996| ) . 



5. INTERPRETING THE TULLY-FISHER RELATION 

We now seek a physical interpretation of the results 
presented in JI] in the context of current models of disk 
galaxy assembly First we discuss various procedures 
for esti mat ing dynamical masses from our rotation curve 
data ((5.1). We then derive estimates of the total bary- 
onic mass (J 5.3). We combine the two estimates to eval- 



uate the relative roles of baryons and dark matter out 
t o the observable radii probed with our deep exposures 
(j |5.4[ ). Although there are considerable uncertainties in 
what follows, our intent at this stage is to illustrate the 
possibilities that will arise when gas masses can be deter- 
mined for samples such as ours so that the total baryonic 
components would be accurately measured and their role 
in the Tully-Fisher relation established. 

5.1. Dynamical Mass Estimates 

The physical basis of our interest in the Tully-Fisher 
relation is that the dynamical mass is strongly corre- 
lated with the luminous and stellar mass components of 
galaxies, and by analyzing empirical constraints, we can 
gain an understanding of the relative assembly histories 
of dark and baryonic matter in galaxies. We thus seek 
to use our data to estimate both the dynamical masses 
(i.e., the total mass, including dark and baryonic) as well 



ture ( 


Pizagno et al. 


2005 


Gnedin et al. 


2007 


et al. 


2010 all low-redshift galaxies) derived ( 



| Williams 
lynamical 

masses trom kinematic data that probe sufficiently far in 
radius to detect the dark halo by revealing a deficit of 
baryons when dynamic masses are compared to stellar 
masses. 

However, our method of using emission line velocities 
to estimate the mass within a given radius, such as T2.2, 
depends sensitively on the assumed shape of the under- 
lying gravitational potential, and hence the distribution 
of mass throughout the disk. For a given ellipsoid po- 
tential, the velocity can be most simply approximated 
as: 

GM(r) 



where £=1 in the case of spherical symmetry. Assuming 
a spherical potential will likely overestimate the disk mass 
unless a spherical dark matter halo is dominant within 
the relevant radius. As such it supplies an effective upper 
limit for a given mass of an ellipsoid calculated from the 
observed circular velocity. Traditionally, dynamical disk 
masses have been calculated with an exponential 'Free- 
man' potential, solved with modified Bessel functions by 
assuming a constant mass- to-light ratio an d an infinitely- 
thin disk of infinite size (Freeman 1970). This ignores 
the presence of the bulge and halo, known t o be impor- 
tant even at the scales considered here (e.g., |Trott et al. 
2010||Dutton et al |[20lTb| > . Therefore, as an alternative 
method of estimating a lower limit, we adopt an oblate 
potential, characterized by a flattening factor q, which is 
the ratio of the scale length normal to the disk over the 
scale l ength of the disk. As shown in |Binney fc Tremaine| 
( |1987[ ) , the velocity for an oblate sphere can then be con- 
sidered as: 



V c {rf w AirGq 



p(m 2 )m 2 dm 
o y/r 2 — m 2 (l — q 2 ) 



(7) 



where m 2 = r 2 + r 2 , and p is the assumed density func- 
tion. The exact shape of the potential will depend on 
the relative contribution of luminous and dark matter, 
as well as on the triaxial shape of each component. Al- 
though halos are believed to be prolate on large scales, 
their shape is less clear at the scales considered here. 
Lensing and dynam ical studies of individual systems 
( Dutton et al.|2011b ) suggest that they may be consider- 
ably rounder. Furthermore, the presence of a bulge gen- 
erally implies that the stellar distribution is significantly 
less flat than that of a pure disk. We adopt q = 0.4 
and an exponential density function as a representative 
maximum oblateness, equivalent to £ s» 0.752 for Eq(6j 
If we were to use a de Vaucouleurs profile (Sersic profile 
where n = 4) instead of an exponential density profile, 
£ would be s» 0.833, resulting in a less than 10% change 
in the dynamical mass calculation. In the following we 
will consider the q = 1 spherical case and the q = 0.4 
oblate case as bracketing the shape of the total potential. 
Additional systematic uncertainties include the effects 
of non-streaming motions, warps and non-gravitational 
forces, as discussed in the well-established literature on 



V c (r) 2 « £- 



(6) 



the interpretati on of local rotation curves (see Binney & 
Tremaine||1987 and references therein). Finally, we con- 
sider possible biases arising from slit spectroscopy and 
its maximum effect propagated to our dynamical mass 
estimates in the following section. 

5.2. Slit-effect correction 

Recent progress with integral field unit (IFU) spec- 
trographs has illustrated some limitations of traditional 
long-slit and multi-slit techniques in determining the in- 
ternal dynamics of intermediate redshift galaxies. This 
slit-effect is similar to beam smearing in radio astronomy, 
where the range of velocities from incoming light are av- 
eraged over the width of the slit, resulting in a broaden- 
ing of Doppler shifted lines and an average reduction of 
the rotational velocity. The magnitude of this effec t has 
been co nsidered in detail b y Kapferer et al. (2006) and 
used by Flores et al. ( |2006 ) to compare IFU -derived ro- 
tational velocities to those determined with a multi-slit 
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instrument. Kapferer et al. (2006) systematically inves- 
tigated the effects oi various slit widths in combination 
with inclination, spatial binning and position angle off- 
sets on measured disk velocities using N-body/SPH sim- 
ulations. We can use these results to consider the effect 
our slit geometry might have in distorting our velocities 
taking into account the galaxy sizes and shapes relative 
to the DEIMOS slits. Kapferer et al. find no systematic 
bias due to binning and posit ion angle offset (beyond 
the correction already made, {3.2). We can, however, 
use Kapferer et al.'s results to calculate an upper-limit 
approximation of the correction to the velocity V2.2 for 
the effect of the slit width relative to scale radius (r s ) 
of each galaxy according to its axis ratio (b/a) and for a 
given inclina tion i. Derived from Figs 9 & 10 of K apferer 
et al. (2006), the correction is: 



Vcorr — Vobs ~t~ Vbe 



1 r s i it 



(b/a) 



sin 2, 



(8) 



20 km s _1 . We find correlations in the 
correction with respect to scale radius and inclination, 
but none with mass, redshift, position angle or observed 
velocity. The correction added to V2.2 ranges from 2 km 
s _1 to 52 km s , with a mean of 20 km s _1 . These 
corrections can be found in our catalog. 

Because of the imprecise nature of these corrections, 
arising from the fact that the Kapferer et al. result as- 
sumes a symmetric Gaussian to the spectral line profile 
(while we fit for two half-Gaussians to account for much 
of the blending between the seeing and dispersion) , the 
Kapferer et al. -based correction remains an upper-limit. 
Thus we did not include them in our precisely measured 
Tully-Fisher relations in §4. However, we will apply them 
to our dynamical mass calculation in order to not bias 
our estimates in a way that may over-estimate the dom- 
inance of the stellar mass compared to the dark matter. 
The difference with and without the slit-effect correction 
can be seen in Fig. [To} Because of the imprecision of 
the analytical formula derived here, we strongly advise 
to not apply such a formula beyond its tested range. 

5.3. Baryonic mass estimates 

In order to examine the redshift-dependent fraction of 
baryonic mass within r^.i we need to obtain an estimate 
of the total b aryo nic mass. In addition to stellar masses, 
discussed in § |2.4| we need to account for the presence of 
gas. 

Accurate gas masses are not yet available for inter- 
mediate redshift galaxies, although CO-derived masses 
have begun to appear for so me systems at z>l with, 
e.g., t he PdBI interferometer ( Tacconi et al.|2010 Daddi 
et al. 2010). Nonetheless the situation will improve sig- 
nificantly through upcoming facilities such as ALMA, 
MeerKAT, and eventually the Square Kilometer Array. 
Although what follows is somewhat speculative, it pro- 
vides a reasonable illustration of what will soon be pos- 
sible. To make progress, we estimated gas masses (M g ) 
for our sample using the local stellar-to-gas mass (M*- 
to-M n ) ratio as a function of recently parameter 
ized bylPeeples & Shankar (2010) based on H I measures 
The H 1 



JNearby Ga 



axy Survey (THINGS) and 



from 

helium-corrected, CO-derived H2 masses from the HERA 
CO-Line EXtragalactic Survey (HERACLES) and the 



Berkeley-Illinois-Marylan d Association Surv ey of Nearby 
Galaxies (BIMA SONG) flLeroy et al.|[2008| . 

According t o " 
ShankarlpOiol): 



the parameterization" By Peeples fc| 



= KfM, 



(9) 



where Kf — 316228, 7=0.57, and M* is measured in 
units of solar masses. 

Until precision gas m asses become available we can- 
not be certain that the Peeples & Shankar (20101 for- 
malism can be applied in this manner at intermediate 
redshift. However, locally- measured gas-to-stellar mass 
ratios will underestimate the gas mass for intermediate- 
redshift galaxies since many stars have subsequently 
formed. To correct for this, we consider the observed 
evolution in the specific star-format ion rate (sSFR) fo r 
blue galaxies to z~2 as measured by Oliver et al. (2010). 
To determine the correction for each galaxy, we integrate 
the best-fit sSFR relation out to the relevant redshift. 



Oliver et al.| (j2010j) find: 

sSFR = X(l + z) a 
where log 10 X/Gyr^ 1 = — 1.36±0.41 and a 



(10) 
-3.4±0.3. 



Madau et al. 


1998; 


Cole et al. 


-40% for a C 


rabrier 

— i rrrr 


q2003p ini- 



2001[ ), estimated to be 
tial mass function from z^ 
is not strictly a closed-box model, where a galaxy sits 
in a huge reservoir of gas and simply converts that gas 
to stars over time. For one, we include the substan- 
tial Chabrier recycling rate consistent with the IMF we 
assume in our stellar mass estimates. The scenario we 
present here is also consistent with roughly equal inflow 
and outflow since z ~ 1 . The fraction of feedback driven 
material that exceeds the escape velocity of the galaxy's 
potential well, and also the typical redshifts at which the 
filaments that feed galaxies evaporate, are still largely 
unknown. The gas mass estimates we add to our stellar 
masses are on average 30% of the total stellar mass esti- 
mates, and only exceed that of the stellar mass in a few of 
the lowest mass objects. So we claim that our gas mass 
estimates are not the main drivers of our conclusions, 
and each galaxy's stellar mass and redshift essentially 
determine the gas mass estimate for this scenario. 

5.4. Comparison of Baryonic and Dynamical Masses 

We finally turn to our comparison of the dynamical 
mass with the stellar and baryonic mas s. The results of 
this comparison are given in Figure 10 The top panels 
show the stellar-to-dynamical mass ratio for a spherical 
and oblate potential. Across our redshift range we find 
our stellar-to-dynamical mass fractions are ~ 0.3 for a 
spherical potential and — 0.4 for an oblate potential, with 
considerable scatter. The dynamical mass estimat es fo r 
the points plotted include the slit-effect correction ({ 5.2 ), 
representing an upper-limit for the dynamical masses, 
and thus a lower limit for the stellar-to-dynam ic m ass 
fraction. For comparison in the plots of Figure [lOj the 
dotted lines are derived from the uncorrected data, and 
represent a lower limit to the dynamical masses, given an 
assumed potential shape (spheroid or oblate in our ex- 
ample). The points without slit-effect corrections are not 
plotted, and just their 3er best-fit contours are plotted for 
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simplicity. In a few cases, the stellar-to-dynamical frac- 
tions are as high as unity, suggesting that baryons play 
a significant role in driving the TF relation. A dominant 
ba ryon fraction is broadly consistent with earli e r results 
bylConselice et al.| J2005D, |Gnedin et al | J2007D, |Dutton 



fc van den ijosch| ( |2U(]gf and [Putton et al.| ( |26lOa| ), es- 



FlG. 10. — Top panels: stellar-to-dynamical mass ratio within our fiducial aperture as a function of redshift, assuming a spherical 
potential (left) and an oblate potential (right). Curves indicate 3<r variations around the best-fit linear relation, with slope b and intercept a. 
The solid lines are derived from dynamical masses including the slit-effect correction, and the dotted lines are derived from the uncorrected 
dynamical masses. The scatter around the best fit relation is also given in each plot. Bottom panels: same as the top panels for baryonic 
(stellar plus gas) mass. 

~ 50% and ~ 70 % respectively for spherical and oblate 
potentials (including slit-effect corrections). It must be 
remembered this baryonic estimate does not include any 
ionized gas. 

Next we explore the radially dependent profile of the 
dynamical mass as calculated from the rotation curve, 
to compare to the baryonic component mass profile. At 
each tenth of a scale radius along the profile, starting 
at 1 scale radius, we compute the dynamical mass given 
the best modeled velocity at that radius in the rotation 
curve, and compare that to the enclosed baryonic mass 
at that radius. The baryonic mass is assumed to follow 
the stellar mass profile, which we approximate to follow 
the distribution of light found in the reddest HST fil- 
ter (F850LP), in an aperture stellar mass estim ate, the 
method of which is described in detail in §2.4| The en- 
closed light at each tenth of a scale radius is compared to 



pecially considering the many uncertainties involved in 
both mass estimates, including the stellar initial mass 
function (IMF). 

The hypothesis that baryonic mass primarily governs 
the slow redshift-dcpcndcnt trends in our observed TF 
relations is supported further when we attempt to add 
estimates for the missing gaseous components to our stel- 
lar masses. We then find that our baryonic mass es- 
timates within r 2 .2 approach those determined dynami- 
cally, with no redshift dependence (bottom panels of Fig- 
ure 10 ). The average baryonic-to-dynamical mass ratio is 
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the total light, stepping out in elliptical apertures from 
the photometry, and the aperture mass is calculated from 
these ratios. 

Although past studies of rotation curves and the TF re- 
lation were often motivated as a means of tracing the in- 
fluence of dark matter on the baryonic assembly of galax- 
ies, ( e.g., |Rubin et al.||1980l [Vogt et al.||1997[ |Pizagno 

20051, it has been suggested 



et al.||20d5| |(Jonselice et al.| 
in local TF studies that the role of dark matter is minimal 



in the optical disk (e.g., Courteau fc R ix 1999; Pahmas & 
|Williams||2000{ |Bell fc de Jong|200l| |Kassin et al.||2006) . 
Figure 1 11 1 shows that, irrespective of the aperture /, de- 
fined in terms of the scale radius (i.e. /=2.2 for r 2 .2), 
out to the furthest observable extent in our disks, the 
mean baryonic-to-dynamical fraction is consistent with 
unity when using an oblate potential for our dynamical 
mass estimates. In contrast, for a spherical potential, 
we find our results permit an equal contribution from 
dark and b aryonic matter by 7-2.2, consistent as we noted 



above wit h Portina ri & Sommer-Larsen (|2007| and Dut- 
ton et al.| ( [2010bp . 



We also show the bary onic-to-dynamic al fraction with 
respect to aperture for a Freeman] (|1970[) potential. The 



unphysical and sharp rise above unity for small radii 
shows that the dynamical mass is likely underestimated 
for the central-most part of the galaxy using this po- 
tential. If used in a maximal disk fit, which would re- 
normalize the sharply rising peak down to unity, the fit 
would likely result in an underestimate of the overall dy- 
namical contribution of the baryons for what is nominally 
the scenario where baryons are maximally contributing 
to the potential. This gives us some indication that the 
Freeman potential, which assumes a constant mass-to- 
light ratio for an infinitely-thin exponential disk, may 
not be the best approximation in maximal/minimal fit 
applications for a disk of a finite thickness and size. 

We thus conclude for the scenario presented here that 
the contribution by baryons to the total mass within the 
radial range probed by observations is between 50% and 
100%. Baryons appear to be the most important compo- 
nent within 2.2 scale radii, and perhaps to larger radii, 
depending on the uncertain conversion from circular ve- 
locity to stellar mass. In addition, such a high fraction of 
baryons is expected to influence significantly the overall 
profile of the dark matter halo at t hese scales through 



gravi tational interactions (see, e.g., Noordermeer et al. 
|2007[ and references therein for a discussion of this topic 
in the local universe). Therefore even if dark matter is 
present in equal amount, it will be tightly coupled to the 
baryonic content. 

6. SUMMARY & DISCUSSION 

Using spectra with extended 6-8 hour integration times 
secured with the DEIMOS instrument on the Keck II 
telescope we have measured 129 rotation curves for 
galaxies selected with zab<22.5 in the redshift range 
0.2<2<1.3 in the two GOODS fields. Using associated 
HST imaging data, we derive modeled velocities cor- 
rected for the effects of inclination and seeing to a fidu- 
cial radius, r 2 .2, corresponding to 2.2 times the disk scale 
length. As 90% of our rotation curves flatten within this 
radius, this gives us a highly reliable set of kinematic 
data which we use to construct the stellar mass (M*) 
and B-band (Mb) Tully-Fisher (TF) relations and their 



M dyn est. from V 
assuming: 

+ Freeman $ 
• Spherical <& 
o Oblate * 



V 



V 



1.0 




1 2 3 4 ! 

Aperture / (r ap /r s ) 

Fig. 11. — Estimates of the mean baryonic-to-dynamical mass 
fraction over all redshifts as a function of relative aperture (/ = 
r/r s ). Shaded regions show the percentage of disks traced to f=2.2 
and 3.2 (r^j); our Keck survey adequately samples the rotational 
velocities to / ~3.5. Different curves relate to different determi- 
nations of the dynamic mass (see text for details) and the tight 
contours show the 1 and 3 times the standard error on the mean. 
Light shading around the results assuming an oblate potential in- 
dicate the aperture-dependent region containing 68% of the indi- 
vidual fractions for our disks. These regions are not shown for the 
other two potentials, but are of similar size to that of the oblate 
potential. 



evolution over the past 8 Gyr. 
results. 



We find the following 



1. We demonstrate a significantly reduced scatter 
around our intermediate redshift TF relations com- 
pared to that seen in earlier studies. The scat- 
ter around the M*-TF relation is typically 0.2 dex 
of M„/M , which is 2-3 times less than earlier 
work, comparable to that determined l ocally and 
less than that achieved bylKassin et al. (20071) who 



introduced an additional dispersion term m order 
to achieve a tight relation. The scatter around 
our Afs-TF relation is typically 0.4-0.7 mag, again 
comparable to that in local relations and a factor of 
2-3 improvement over earlier work at intermediate 
redshift. In addition to demonstrating the validity 
of our observational approach and our improved 
modeling techniques, our results clearly show the 
TF relation is well-established at redshift z ~1. 

2. We find a modest but statistically-insignificant evo- 
lution in the M*-TF relation with redshift corre- 
sponding to a growth in stellar mass at fixed veloc- 
ity of A M* = 0.04 ±0.07 dex with cosmic time 
from a redshift (z)~1.0 to 0.3. This is consis- 
tent with the results of recent numerical and semi- 



analytic simulations dPortinari & Sommer-Larsen 



2007j |Somerville et al!pU°U8| |Dutton et al.||2011a| ). 



3. More pronounced evolution is seen in the M^-TF 
relation corresponding to a decline in luminosity of 
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0.85 ± 0.28 mag from (z)~ 1.0 to 0.3, again con- 
sistent with p redictions from Portinari & Sommer- 



Larsen 



( |2007 [). The origin of this evolution can 
be understood by examining the mass-dependent 
correlation between M„ and Mb as a function of 
redshift. Together with the larger scatter seen in 
the Mb-TF relation than in that based on stellar 
mass, this demonstrates that the M*-TF relation 
is the more fundamental descriptor of disk galaxy 
assembly. 

4. We finally use our data to illustrate the poten- 
tial of comparing dynamical and baryonic masses 
to better understand the contributions that they 
matter make in the M*-TF relation. Although we 
are uniquely placed to explore our extended rota- 
tion curves and carefully-derived stellar masses, our 
gas mass estimates are clearly approximate at this 
stage. We estimate dynamical masses using both 
spherical and oblate potentials which we expect 
bracket the likely values. We compute baryonic 
masses by estimating the additional mass in cold 
gas. We find that baryons may contribute between 
50-100% of the total mass within the radii at which 
we can confidently observe dynamics of the gas in 
our disk sample. Such a high fraction of baryons 
influence significantly the overall profile of the dark 



matter halo. Even if dark matter is present in equal 
amount, it will be tightly coupled to the baryonic 
content. 
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